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Abstract 

It is known experimentally that stable few-body clusters containing 
negatively-charged electrons (e) and positively-charged holes (h) can exist 
in low-dimensional semiconductor nanostructures. In addition to the famil- 
iar exciton (e+h), three-body 'charged excitons' (2e+h and 2h+e) have also 
been observed. Much less is known about the properties of such charged 
excitons since three-body problems are generally very difficult to solve, even 
numerically. Here we introduce a simple model, which can be considered as an 
extended Calogero model, to calculate analytically the energy spectra for both 
a charged exciton and a neutral exciton in a one-dimensional nanostructure, 
such as a finite-length quantum wire. Apart from its physical motivation, 
the model is of mathematical interest in that it can be related to the Heun 
(or Heine) equation and, as shown explicitly, highly accurate, closed form 
solutions can be obtained. 
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I. Introduction 

The optical properties of low-dimensional semiconductor structures, called nanostruc- 
tures, have attracted much attention in the past few years, both experimentally and theo- 
retically One of the most interesting questions concerns the properties of 'excitons' in such 
low-dimensional structures. An exciton (X) is a neutral, two-body complex formed by the 
attractive force between a negatively- charged electron (e) in the semiconductor conduction 
band and a positively-charged hole (h) in the valence band. An exciton therefore appears 
to be somewhat analogous to a hydrogen atom. There are, however, two important dis- 
tinctions. First, the hole and electron masses are typically of the same order of magnitude 
and, second, the low dimensionality of the nanostructure can restrict the electron and hole 
motion to such an extent that the exciton must be treated as either two- or one- dimensional. 

Recent observations of anomalies in the optical spectra of quantum wells (i.e. two- 
dimensional nanostructures) have been attributed to the formation of negatively- charged 
excitonsR Such complexes can arise when an exciton is created in the presence of a low 
concentration of free electrons; it may then be energetically favorable for the exciton to 
capture one electron to form a negatively-charged exciton, i.e. X + e — > X . In addition 
to two-dimensional nanostructures, it is interesting to consider the possibility of X~ for- 
mation in one- dimensional nanostructures such as a quantum wire. The consideration of 
such low- dimensional systems is particularly important since the exciton binding energies 
in a quantum wire are higher than those in the quantum well, owing to the reduced dimen- 
sionality; this increased exciton binding energy is thought to underlie the recently observed 
exciton lasing in a quantum wire devicei. 

This paper uses a simple model to investigate the properties of three-body complexes such 
as charged excitons. In particular, we provide closed-form expressions for the energy spectra 
of both a charged exciton X~ and a neutral exciton X in a finite-length, one-dimensional 
quantum wire. Even in one dimension, the three-body problem (i.e. two electrons (e) and 
one hole (h)) with a Coulomb interaction would necessitate a computationally-intensive 
numerical solution. This is not our goal; instead we wish to demonstrate that analytically- 
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solvable models can be introduced to identify trends in the X - and X energy spectra as 
a function of device parameters. We therefore sacrifice quantitative accuracy concerning a 
particular device in favour of a broader understanding of qualitative behaviour. 

Our model considers an inverse-square interaction potential between particles (e-h is 
attractive and e-e is repulsive) together with a parabolic confinement potential of arbitrary 
strength along the wire. A non-Coulombic form for the interaction is in fact not unrealistic 
in nanostructures due to the presence of image charges in neighbouring gates and electrodes 
(see, for example, Ref. 3). As will be discussed, the energy spectrum of the neutral exciton 
X shows the same qualitative behaviour with both 4j (i.e. inverse-square) and (i.e. 
Coulomb) interactions. For X - , the three-body Schrodinger equation is shown to reduce 
to the Heun equation. The complete energy spectrum is found in the regime of physical 
interest. The analysis suggests that the X complex can have an enhanced stability as 
compared to X. 



II. Neutral exciton X 

Our model Hamiltonian for the neutral exciton X (i.e. e-h pair) in a finite-length one- 
dimensional quantum wire is given by 

„ _ _ h 2 ( d 2 d 2 \ 1 , 2f 2 , 2 ._ g eh h 2 1 m 
H ~ 2m*{dxl + dxl) + 2 mUJ ° [Xe + Xh) 2m*{x e -x h y Uj 

where x e and x h are the electron and hole coordinates, the dimensionless parameter q eh 
characterizes the electron-hole interaction strength, and m* is the effective mass of the 
electrons and holes (assumed identical). The parabolic confinement potential has arbitrary 
strength and is assumed to be the same for both the electron and hole; the confinement 
parameter ujq can be chosen so as to mimic the effect of a wire of finite length L since 
L 2 ~ ft(m*u;o) _1 - If q eh > | in Eq. (1) then it is trivial to show that the two particles collapse 
toward each other for any finite value of E; this can be seen by examining the behavior of 
the wave function near x e ~ Xh- To avoid this we assume that < q e h < \- In terms of 
the center-of-mass coordinate X = \{x e + Xh) and the relative coordinate x = x e — Xh, the 
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Hamiltonian becomes H = H c . m ,(X) + H re i(x), where H c , m , is the Hamiltonian for a single 
particle in a one- dimensional harmonic potential and 

^ d 2 1 2 2 q eh h 2 1 

with a reduced mass /i = The eigenvalues of i? c . m . are the one- dimensional harmonic 
oscillator levels. The eigenvalues of H re i are! 
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E rel (n;q eh ) = 7kJ Q (2n+ - - A) ; n = 0, 1, 2, . . . (3) 



where 



A = i-i(i-2 9efc )i . (4) 

As the e-h interaction is reduced (i.e. g e h — > 0) Eq. (|3|) becomes 
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£(n; 0) = hw (2n + -) ;n = 0,1,2,... (5) 

i.e. the odd energy levels of a harmonic oscillator!. We shall refer to the quantity A as the 
electron-hole energy-shift for reasons which are clear by comparing Eq. (Q) and Eq. (||). 

The Hamiltonian for the neutral exciton with a Coulomb interaction, Xc, is the same 
as in Eq. ([I]), except that the last term is replaced by the interaction potential r^. We 
emphasise, however, that in the presence of image charges in neighbouring gates etci, a 
Coulomb interaction will not necessarily be more realistic than an inverse-square interaction. 
To date, the one- dimensional Schrodinger equation for Xc has not been completely solved 
analytically. The Hamiltonian for has the form 

rrC ^ & 1 2 2 ^ 1 

where e e h is the dielectric constant of the system. 

Introducing the dimensionless coordinate y, where x = by with b = the eigenvalue 

problem of Eq. @ reduces to: —i/}"(y) + {y 2 - ^^ipiv) = ^E re i/huJoip(y). Equation (H) 



reduces to: —ip"{y) + (y 2 — f^v/paTipV'Cz/) = ^rei/^^iv)- Typical values for the device 



parameters are as follows: e e h = 12, Tiloq = O.OleV and m*/m = 0.07. Using these values, 
we have performed numerical calculations which show that the low-lying energy spectra for 
Xc and X can indeed be quantitatively similar, provided an appropriate value of the free 
parameter q e h is chosen. In the q e h — > limit Xc yields an identical energy spectrum to that 
given in Eq. (Q). This indicates that the inverse-square interaction is 'as realistic' as the 
bare Coulomb interaction in the context of one- dimensional models. To understand why this 
should be so, one can consider trying to solve the X and Xc problems using the complete 
basis set of the one- dimensional oscillator eigenstates \n) in the interval — oo < x < oo, i.e. 

(x\n) = A n E n (y)exp(-y 2 /2) 

where y = ^Jmuo/hx; H ra is a Hermite polynomial and A n is the appropriate normalisation 
constant. In the complete basis of kets \n) we need to evaluate matrix elements of the type 
(m\H\n), where H is either the Hamiltonian of X or Xc- This matrix element includes terms 
like 

/oo 
R n (y)R m (y)V(y)exp(-y 2 )dy 
-oo 

where the potential V(y) is either proportional to ^ or ^. It is clear that the integral 
diverges when m and n are both even. As discussed in Ref. 6, this integral is only finite if 
we truncate our Hilbert space so that m and n are odd; this finding is hence consistent with 
the result of only odd energy levels in Eq. In this sense, both the A and ^ potentials 
are 'non-penetrable'@. Physically, this means that the electrons and hole in the X and Xc 
complexes cannot interchange their particle ordering along the wire; the exchange energy is 
therefore zero. The configuration with x > is totally separate from the configuration with 
x < 0. Each energy level of X and Xc is therefore doubly degenerate when the complexes 
are defined in the full interval — oo < x < oo. Equation (|3|) represents the complete energy 
spectrum of the neutral exciton X with inverse-square interaction. 

III. Charged exciton X~ 
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Our proposed Hamiltonian for the X~ complex may be considered as a generalisation 
of the Hamiltonian discussed by Calogeroi. Calogero considered the three-body problem 
with a harmonic oscillator potential and inverse-square pair potentials for the case of three 
identical particles. Our model Hamiltonian H for two electrons and one hole has the same 
form, but we allow the strength and sign of the interaction between the particles to be 
different. In particular, 

i=el,e2,n \ t / 

qehh 2 1 | q ee h 2 1 q e hh 2 1 ,^ 

2m* (x e i — Xh) 2 2m* (x e \ — x e 2) 2 2m* (xh — x e 2) 2 

where x e i, x e 2 and Xh are the coordinates of the two electrons and hole, q ee is the electron- 
electron interaction parameter (q ee > 0); as in Sec. II we restrict the electron-hole in- 
teraction parameter to < q G h < §. The three-body problem in Eq. (7) is separable^. 
The separation involves two coordinate transformations; first a Jacobi transformation: 
X = \{x e \ + x e 2 + Xh), x = 2~ x l 2 [x e \ — Xh) and y = 6~ 1 ^ 2 (x e i + Xh — 2x e2 ) which enables us 
to rewrite H = H cm (X) + H re i(x,y). Second, we write H re[ in terms of polar coordinates 
x = r sin(<f) — 7r/3) and y = r cos(<f) — vr/3): 

TT , „ ft 2 ( d 2 I d\ I 22 h 2 M 

H rel (r, 0) = -— ^— + - - j + -m u r + —-, (8) 

where M is the following operator 
d 2 1 " 



M tkfi ' 2 



q e h qee q e h 



(9) 



sin 2 (cf) — n / 3) sin 2 ((J)) sin 2 (<f) + tt / '3) 

Writing an eigenfunction of Eq. (8) as ip(r,(f)) = i?(r)$(0), the three-body problem has 
therefore been reduced to finding the solutions of the following two ordinary, second-order 
differential equations 

H rel R(r) = E rel R(r) (10) 

and 



M$ z (0) = (11) 
The eigenvalue problem in Eq. (10) is solved in Ref. 4 and the eigenvalues are given by 

E rel = %uj {2n + 1 + bi) ;n = 0,1,2,... (12) 



leaving the nontrivial problem of solving Eq. (pTTj) . For a given particle configuration the 
angle <fi is limited to a certain interval. The ordering of the three particles (i.e. 2e + 
lh) is therefore determined by <fi. As seen earlier for X, the 'non-penetrable' interaction 
potentials prevent particle interchange and therefore make it necessary to treat different 
particle configurations separately. In the interval G]0;7r[ the relationship between and 
the particle configuration is given by 

4> G]0; 7r/3[ : x e2 < x el < x h 

<\> g]tt/3; 2tt/3[ : x e2 < x h < x el (13) 
4> €]2tt/3; 7r[ : x h < x e2 < x el . 

Using the trigonometric identity sin3<f) = 3sin<j) — 4sin 3 <f) we can rewrite the wave function 
in Eq. (11) as 

$(x) = (x- l) A - /2 (x - ^ h y(x) (14) 

where x = cos 2 ((p), A eh = 1 - A = 1/2 + 1/2(1 - 2q eh ) l l 2 and A ee = 1/2 + 1/2(1 + 2g ee ) 1 / 2 . 
Hence we have reduced the problem of solving Eq. (11) to solving the Heun differential 
equation! 

d 2 V (l 5 e \ dy \ af3x - q =0 (15) 

dx 2 \x x — 1 x — l/AJdx x(x — l)(x — 1/4) 

with coefficients 

7 = 1/2 a = A ee /2 + A eh ,-l/2 6 ; 

S = A ee + 1/2 /3 = A ee /2 + A e ^ + 1/2 k (16) 

e = 2A e , q = l/2(A eh ) 2 + (A ee /4) 2 - (6 z /4) 2 

and with the five parameters satisfying 



a + f3 -7 - 5 - e+1 = . 



(17) 



q is the so-called accessory parameter. An introduction to the general features of Heun's 



Heun equation called the Heine equation^. To our understanding, Heine's equation is far 
less well-known and we will therefore continue referring to Eq. (|T5| ) as the Heun equation. 
Solutions to Eq. (|l^) which are of particular interest are those which are analytic in some 
domain enclosing two singularities. Such solutions are called 'Heun functions' and are de- 
noted by Hi, using the notation of Ref.i. We employ this notation in order to distinguish 
Heun functions Hi from Heun polynomials Hp, which are analytic in an interval containing 
three singularities. Below, Heun functions which are analytic in the interval x G [0; 1/4] or 
G [7r/3;27r/3] will be studied in more detail. An eigenfunction of Eq.(]T5|) defined in the 
region G [tt/3; 2tt/3] can be written, using cos 2 (ft = 1 — sin 2 (ft, as 



Heun functions may be found by the power-series method or by the method of hyperge- 
ometric function seriesi. However for both methods the coefficients in the series have to 
satisfy a three-term recursion relation. It might be possible to use the method of augmentet 
convergence to extract Heun functions from such relations, but in most cases only a numer- 
ical procedure is possible^. Although we cannot in general solve the eigenvalue problem of 
Eq. (15) analytically, we will now obtain a set of highly accurate, approximate solutions for 
the region G [7r/3; 27r/3]; these solutions are essentially exact for a large interval of the 
ratio q ee /(leh, including the range of physical interest. Employing the following trigonometric 
formula 



equation is given in Ref.i. In addition, a recent bibliography containing 300 classified entries 
related to Heun's equation are listed in Ref.0. The relation in Eq. (|T7| ) ensures that the 
four singularities of Eq. (|1|) stay regular. One of these regular singularities is, however, an 
elementary singularity (since 7 = 1/2) which implies that Eq. ([15]) is a special case of the 



$;(0) = (sm0) A -(sm 2 0-3/4) Aeh iff(l/4,g;a,/3,7,5; cos 2 (j)) 



(18) 



/ 



+ 
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f 




9/ 



(19) 



sin 2 (0 — 7r/3) 



sin 2 {4>) 



sin 2 (0 + 7r/3) 



sin 2 (30) 
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where / and g are arbitrary functional^, the operator M in Eq. ( |TTD reduces to 

ds 2 2 



(20) 



sin 2 (f) sin 2 3(J) 

where k = q ee /qeh- Consider the angular dependence of the total interaction potential of the 
X complex: V((f)) = (1 + k) / sin 2 (ft — 9/sm 2 30. Here V((f>) is periodic in 0, V((p) = V((p + n), 
hence the physics of the three-body problem is contained within a 0-interval of n, e.g. 
e]0; 7r[. Figure 1 shows V(0) for the case k — 1. The asymptotic behavior of V(4>) is easily 
understood: at = the repulsion between the two electrons causes a positive singularity, 
at = 7r/3 the attraction between the electron at x e \ and the hole at Xh causes a negative 
singularity, etc. For the particle configuration where the hole is between the two electrons 
(in Fig. 1: G]7r/3; 27r/3[ ), the potential does not have a repulsive divergence and therefore 
corresponds to the most stable three-body configuration. A very good approximation to 
V(4>) in this interval (0 e]7r/3; 27r/3[) for moderate values of k is 



^ ) = 1 + K -sJ^ • (21) 



This is illustrated in Fig. 2. Evidently the approximation becomes more exact for smaller 
k and, as will become clearer later, exact in the limit q ee — > —q e h- We shall restrict k to the 
interval k G [0;20[; V app (<j)) is now an excellent approximation to V(<f)). The quantity k is 
the ratio of the strength of the electron-electron interaction to the strength of the electron- 
hole interaction, and for all practical electronic devices we expect this ratio to be less than 
20. Using this approximation, Eq. ( pTTj) becomes M app $;(0) = b 2 §i(4>) where M app = 
—d 2 /d(p 2 + ^YV app . The task of solving this differential equation can now be transformed 
into the problem of solving a hypergeometric equation. The solutions are 

$i(0) = (sm30) Ae,l 2 iq(a, b; c; cos 2 30) (22) 

where 

c = 1/2 , a = l/2(A eh -b[) , b= l/2(A eh + b[) (23) 
and b' 2 = l/9bf — 1/18(9^ + q ee ). The exact eigenvalues are found to be 



k = [9(1 + 1 - A) 2 + ~(q ee + q^)} 1 ' 2 ; 1 = 0,1,.... (24) 

Combining Eq. (^3) and Eq. ( |T2"D gives the complete energy spectrum of the X~ complex for 
the particle configuration where the hole is placed between the two electrons and k G [0; 20 [; 
in particular, 

E re i(n, I- q eh , q ee ) = hu hn + 1 + [9(1 + 1 - A) 2 + l -(q ee + q eh )]^ . (25) 

In the Calogero- model limit where q ee — > —q e h, the energy spectrum in Eq. (l2~5f) reduces 
to the energy spectrum found by Calogero for three identical particles^. By comparing Eq. 
(f22|) and Eq. (|1|) the approximation corresponding to Eq. ([H]) can be expressed as 

$j(0) = (sin(f>) Aee (sin 2 <p- 3/4) Aeh iJf(l/4, g; a, (3, 7, 5; cos 2 <p) 

~ (sin30) Aek 2 Fi(a, 6; c; cos 2 30) (26) 

In the Calogero-model limit, Eq.pB]) becomes exact and Hf = (— 4) Aeh 2-Fi- 

There exist two distinct particle configurations in which the hole is positioned between 
the two electrons. These two configurations represent different physically-accessible systems, 
because of the 'non-penetrable' property of the interaction potential as discussed in Sec. 
II. If X~ is defined for the full range of these two particle configurations, then each level 
in Eq. ( ^5f) is doubly-degenerate. By forming suitable linear combinations, the subspace 
corresponding to each eigenvalue may be spanned by an antisymmetric and symmetric wave 
function with respect to interchange of the two electrons. 

IV. Comparison of energies of X and X 

We have derived the energy spectra for X and X - . One possibility might be to compare 
their relative stabilities by turning off the confinement potential and calculating the binding 
energies of X and X~. However, a finite confinement potential is needed within the present 
model in order to produce discrete energy levels for the complexes; when the confinement 
interaction is turned off (uq — > 0) a continuous spectrum0 is obtained for both X and X~. 
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This suggests than X and X are not exciton complexes in the usual sense, since their 
existence depends on the presence of a confinement potential. Keeping the confinement 
potential finite, the ground state energies of X and X~ in the non-interacting limit (q ee , q e h — > 
0) are given by {E re i + E CmTn ) which yields 2huo and 9/2huo respectively; these energies are 
identical to the ground state energies for two and three spinless fermions in a harmonic 
potential well. We emphasize that this particular non-interacting limit is reached as a 
consequence of the model being strictly one-dimensional and containing a singular potential. 

Keeping the confinement potential finite, we now investigate the changes in energy of 
X and X - as the two-body interaction is varied. We introduce a quantity which we refer 
to as the 'interaction energy' E int , defined as the energy obtained by subtracting the total 
energy with vanishingly small interactions (i.e. q ee , q e h — > 0) from the total energy with 
finite interactions. For X and X in their respective ground-states, we obtain 

Ef nt {q eh ) = E rel (0; q eh ) - E rel (0; 0) = -tuv A (27) 

and 

Ef; t {q e e, Qeh) = E rel (0, 0; q ee , q eh ) ~ E rel (0, 0; 0, 0) = hcu ([9(1 - A) 2 + l -{q ee + q eh )] 1 ' 2 - 3) . 

(28) 

As expected, both Ef nt and E^ t become increasingly negative with increasing q e h, i.e. the 
ground-state energy decreases as the electron-hole interaction q eh increases. In addition 
E^ t increases as the electron-electron interaction q ee increases. We therefore interpret the 
quantities E^ t and Ef nt as indicative of the binding strength of X and X - respectively. 
As a consequence, the relative stability of X~ and X is then effectively represented by the 
'relative binding strength' 

AE int (q ee: q eh ) = Ef~ - Ef nt = hw ([9(1 - A) 2 + l -(q ee + q eh )} 1 ' 2 - 3 + a) . (29) 

If AE int < 0, this would suggest that X~ is more strongly bound than X; if AE int > 0, the 
reverse is true. The cross-over occurs when AE int = 0, i.e. when 
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q ee = 8A + 7q eh . (30) 

We label the k value for which Eq. ( p0|) is satisfied to be K eq . When k < K eq for a given 
q e h, then AEi nt < 0; when k > K eq the reverse is true. This is illustrated in the inset of 
Fig. 3 where AE is plotted as a function of k for the case q e h = 0.25. The main part of 
Fig. 3 shows function of q eh . Two features are interesting. First, n eq stays within 

a reasonably narrow interval 11 < n eq < 15 for all values of q e h within the model. In the 
range < q e h < 0.4, K eq increases only slightly on an absolute scale. Second, the minimum 
value for n eq is given by K eq = 11. Therefore if k is smaller than 11, as is the case for typical 
experimental devices, this would suggest that X~ is more strongly bound than X. 

V. Conclusions 

In summary, we studied a new model for a class of three-body problems; highly accurate, 
closed form solutions for this model were obtained. The model was used to study the exciton 
complex X~ relative to the neutral exciton X. The analysis suggest that X~ might be more 
strongly bound than X for typical one-dimensional devices. 

Finally we note that the present results for X~ apply equally well to X + (i.e. one electron 
plus two holes). Future work will use the same model to examine more exotic excitonic 
complexes containing n e electrons and holes, i.e. X An ~ where An = n e — rih- 
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Figure Captions 

Figure 1. 

The potential function V{4>) for the case k = 1. Also shown schematically are the 
corresponding electron-hole configurations corresponding to the various 0-regions. 

Figure 2. 

Comparison of the approximate potential V app ((f)) and the exact potential V(<f>) for various 
values of k. For the case k— 1, V(<f>) and V app (4>) are essentially identical. 

Figure 3. 

The ratio of the electron-electron interaction to the electron-hole interaction at the sta- 
bility cross-over point, function of the electron- hole interaction q eh . The inset shows 
AE int as a function of k for the case q e h = 0.25. 
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